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With our recently proposed effective Hamiltonian via Monte Carlo, we are able to compute low energy physics 
of quantum systems. The advantage is that we can obtain not only the spectrum of ground and excited states, 
but also wave functions. The previous work has shown the success of this method in (l+l)-dimensional quantum 
mechanical systems. In this work we apply it to higher dimensional systems. 



1. Introduction 

Monte Carlo (MC) method with importance 
sampling is an excellent technique to calculate 
high dimensional integrals, and has successfully 
been applied to Lagrangian lattice gauge theory. 
Recently, we proposed an algorithm to compute 
the MC Hamiltonian The advantage, in com- 
parison with the standard Lagrangian MC ap- 
proach, is that one can obtain not only the spec- 
trum of ground and excited states, but also the 
wave functions. The method has been tested 
in quantum mechanics (QM) in 1+1 dimensions, 
where the exact results were reproduced with high 
precision jl] M. In this paper, we apply this algo- 
rithm to QM in 2+1 dimensions. The motivation 
is as follows, (a) A number of algorithms for solv- 
ing the Schrodinger equation in lower dimensions 
(e.g. Runge-Kutta) do not work in 2+1 dimen- 
sions. Our method works in 2+1 dimensions, it is 
likely to work also in higher dimensions, (b) The 
ultimate goal of the algorithm is to study struc- 
ture function, the S matrix, QCD at finite density 
as well as systems with many degrees of freedom. 
QM in 2+1 dimensions is a simple example. 



2. Algorithm 

Let's review briefly the basic ideas |Q. Ac- 
cording to Feynman's path integral approach to 
QM, the (imaginary time) transition amplitude 
between an initial state at position Xi, and time 
ti, and final state at Xf, tf is related to the Hamil- 
tonian H by 

< x f ,tf\xi,ti >=< x f \e- H(t f- ti)/n \xi > 



= Y J <Xf\Ev>e- E » T/h <E„\ Xl>1 (1) 

where T = tf — ti. The starting point of our 
method, as described in detail in |jj is to con- 
struct an effective Hamiltonian H e f / (finite NxN 
matrix) by 

< x f \e- H ^-^/ n \xi >«< xt\e- H 'f' T / H \xi > 

N 

= J2< *f\K ff > e~ E l" T ' h < E°Jf\ Xi > . (2) 

H e ff can be found by MC simulation using the 
following procedure: 
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(a) Discretize the continuous time. 

(b) Generate configurations [x] obeying the Boltz- 
mann distribution 



P(x) 



exp(-S[x]/h) 



f[dx] exp(-S[x]/h) ' 
(c) Calculate the transition matrix elements 

M fi =< x f \e- H °" T / h \xi > 



(3) 



(4) 



between N discrete Xi points and N Xf points. 

Note that the matrix M is symmetric. 

(d) Diagonalize M by a unitary transformation 



M 



U^DU, 



(5) 



where D = diag{e~ E i" T / h , e~^s tT / n ). Steps 
(a) and (b) are the same as the standard MC 
method. Step (c) is the essential ingredient of 
our method, from which we can construct H e ff, 
and obtain the eigenvalues El/f and wavefunc- 
tion \E%>f > through step (d). Once the spec- 
trum and wave functions are available, all physi- 
cal information can be retrieved. Since the theory 
described by H is now approximated by a theory 
described by a hnite matrix H e f f , the physics of 
H and H e f t might be quite different at high en- 
ergy. Therefore we expect that we can only repro- 
duce the low energy physics of the system. This 
is good enough for our purpose. 

3. Results 

We consider the following examples in D = 2: 

3.1. Uncoupled harmonic oscillator 

The Euclidean action is given by 



5' 



dt 



-(x 2 +y 2 ) + —(x 2 



y 



(6) 



The spectrum is degenerate and exactly known: 
E ni ,n 2 = ujh(ni + n 2 + 1), rci.,712 = 0, 1, 2, . . . (7) 

3.2. Coupled harmonic oscillator 

The Euclidean action is given by 



S= [ dt 
Jo 



■(x 2 +y 2 ) + Xxy 



The spectrum is also exactly known: 
E nu n 2 = Wi^(m + -) + L0 2 ti(n2 + -), 



711, 7l 2 = 0, 1,2, 



(9) 



where u)\ — 



— and u>2 



uj 2 - —. 



3.3. Summary 

Tabs. 1 and 2 compare the spectrum from H e ff 
with the exact results. Figs. 1-4 show the first 
two wave functions. They are in very good agree- 
ment with the exact ones, at least in the low en- 
ergy domain. In Ref. pj, more results for wave 
functions of higher execited states, and thermo- 
dynamical quantities such as the average energy 
and specific heat, were reported. However, the 
CPU time is much longer than in 1-dimensional 
systems. We need to find a better (stochastic) 
way to select the basis in Hilbert space in step 
(c) of Sect. 2, before extending the method to 
quantum field theory. We are curious to see if our 
method can resolve the long standing problem of 
QCD at finite density. It is worth mentioning the 
recent analytical p] and numerical [pj efforts in 
Hamiltonian formulation. 
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Table 1 

Spectrum of uncoupled hamonic oscilator. 



Tabic 2 

Spectrum of coupled hamonic oscilator. 



m,n 2 


pMC 


Tphlxact 
n, n 9 , 




ni,n 2 


pMC 

n i n? 


TphJxact 


0, 


0.9975 ± 0.0258 


1.0000 




0, 


0.9640 ± 0.0245 


0.9659 


0, 1 


1.9902 ± 0.0420 


2.0000 




0, 1 


1.6693 ± 0.0335 


1.6730 


1,0 


2.0022 ± 0.0468 


2.0000 




1,0 


2.1843 ± 0.0535 


2.1907 


0, 2 


2.9499 ± 0.0393 


3.0000 




0, 2 


2.3650 ± 0.0363 


2.3801 


2,0 


2.9815 ± 0.0550 


3.0000 




1, 1 


2.8818 ± 0.0553 


2.8978 


1, 1 


3.0014 ± 0.0546 


3.0000 




0, 3 


3.0883 ± 0.0447 


3.0872 


n-i + n 2 = 3 


3.9718 ± 0.0496 


4.0000 




0, 3 


3.0883 ± 0.0447 


3.0872 


n 1 +n 2 =3 


3.9853 ± 0.0626 


4.0000 




2,0 


3.4049 ± 0.0532 


3.4154 


n-i + n 2 = 3 


4.0547 ± 0.0611 


4.0000 




1, 2 


3.6275 ± 0.0526 


3.6049 


n\ + n 2 = 3 


4.0713 ± 0.0851 


4.0000 




0,4 


3.7414 ± 0.0327 


3.7944 
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Figure 1. ip niri2 (x,y) vs. x at y = — 1 for uncou- 
pled harmonic oscillator with n\= n 2 = 0. 
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Figure 3. 4> ni n 2 { x i v) vs - x at y = —1 for coupled 
harmonic oscillator with n\ = n 2 = 0. 
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Figure 2. 4Jnm 2 { x -,y) vs - x dX y = —1 for uncou- 
pled harmonic oscillator with n\ = 0,n 2 = 1. 



Fi gure 4. ipnm-zfai y) x c d,t y — 1 for coupled 
harmonic oscillator with m = 0, ri2 = 1. 



